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Abstract 

The Hamiltonian of the spinless relativistic Coulomb problem combines the standard 
Coulomb interaction potential with the square-root operator of relativistic kinematics. 
This Hamiltonian is known to be bounded from below up to some well-defined critical 
coupling constant. At this critical coupling constant, however, the differences between 
all analytically obtainable upper bounds on the corresponding energy eigenvalues and 
their numerically determined (approximate) values take their maxima. In view of this, 
an analytical derivation of (not so bad) upper bounds on the lowest-lying energy levels 
at the critical coupling constant is presented. 
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The "spinless relativistic Coulomb problem" to be investigated here 
is defined by a self-adjoint Hamiltonian H composed of the square-root 
operator of the relativistic expression for the free (so-called "kinetic") 
energy 

T = V^p2 + m2 (1) 

of some particle of mass m and momentum p as well as the (spherically 
symmetric) Coulomb interaction potential 

Vc(r) = --, «>0, (2) 
r 

depending merely on the radial coordinate r = |x| and involving some 
dimensionless coupling constant a: 

H = Jp^ + m^--. (3) 

r 

The serious investigation of the spectrum of the operator (|^) started 
with Ref. Ilj], where it was shown that, roughly speaking, this operator 
is bounded from below precisely up to a critical value ac of the coupling 
constant a, namely, up to 

2 

ac = - . 

TT 

Unfortunately, until now this operator withstood all efforts to localize 
its (discrete) spectrum analytically. (Some brief account of the history 
of these attempts may be found, for instance, in Ref. []2j].) In fact, from 
the point of view of analytical statements about the eigenvalues of this 
operator, at present one has to be content with sets of upper bounds on 
the latter By comparing with numerical computations, however, 

one realizes that, rather generally, the quality of all these upper bounds 
decreases with increasing coupling constant a. In view of the above, in 
this note we are interested in the analytic calculation of upper bounds 
on the eigenvalues of the Hamiltonian H given in (^ for large values of 
the coupling constant a, in particular, at the critical coupling constant 

The theoretical basis as well as the primary tool for the derivation of 
rigorous upper bounds on the eigenvalues of some self-adjoint operator 
is, beyond doubt, the so-called "min-max principle" [^. An immediate 
consequence of this min-max principle is the Rayleigh-Ritz technique: 
Let H be some semibounded self-adjoint operator [which, according to 
a comprehensive spectral analysis presented in Ref. [jlj] , obviously holds 

1 See also Ref. g. 
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for our (semi-)relativistic Hamiltonian, Eq. (^]. Let Ek, /c = 0, 1, 2, . . ., 
denote the eigenvalues of H, ordered according to Eq < Ei < E2 < . . .. 
Let Dd be some d-dimensional subspace of the domain of H and let E^, 
/c = 0, 1, . . . , d — 1, denote all d eigenvalues of this operator H restricted 
to the space D^, ordered according to Eq < Ei < . . . < E^-i- Then the 
kth eigenvalue (counting multiplicity) of iif, E^^ satisfies the inequality 

Ek<Ek, A; = 0,l,...,d-1 . 

Now, let us assume that this (i-dimensional subspace is spanned 
by a set of d linearly independent basis vectors 1?/^^), /c = 0, 1, . . . , d — 1. 
Then the set of eigenvalues E may immediately be determined as the d 
roots of the characteristic equation 

det {mH\^,) - E{^,\^,)) = , z, J = 0, 1, . . . , - 1 , (4) 

as becomes clear from the expansion of any eigenvector of the operator 
H in terms of the basis vectors \il)k)i = 0, 1, . . . , (i — 1, of the subspace 
Dd. 

The Coulombic Hamiltonian (3) involves only a single dimensional 
parameter, namely, the particle mass m. Consequently, for a vanishing 
particle mass all energy eigenvalues Ek must vanish too, as is also seen 
by application of the "relativistic virial theorem" proven in Refs. [0, ^. 
Accordingly, in order to remain on the safe side, we consider here only 
the special case of a nonvanishing particle mass m: 

m > . 

Our particular choice for the basis vectors \iljk) to be adopted here is 
defined by trial functions il) given in configuration-space representation 
(r = |x|) by9 

^^(r) = iVfcr^+^-iexp(-mr) , > , 
with the normalization factor 

(2^)2^+2/3+1 

47rr(2/c + 2/3 + 1) ' 

^ A proper orthogonalization of this basis would yield basis functions which involve 
the (generalized) Laguerre polynomials. In some respect, this orthogonal basis appears 
to be more convenient for a numerical (which means, not entirely analytical) treatment 
of the present problem p. 
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or in momentum-space representation (]7 = |p|) by 



~ _ ~ sin[(A: + arctan(j^/m), 
W)-Nk ^(^2 + ^2)(fc+/5+i)/2 ' 

with the normalization factor 
Nk = 



^ r(/c + /3 + i 



^27r2r(2A; + 2/3 + l) 

both of which are, of course, related simply by Fourier transformation. 
Whereas k indicates a positive integer, /c = 0, 1, 2, . . ., the parameter P 
is introduced to allow, for some given value of the coupling constant a, 
of a total cancellation of the divergent contributions to the expectation 
values of kinetic energy T and interaction potential Vc{r): (3 = (3{a). 
More precisely: in order to provide for that cancellation, the parameter 
(3 must be adjusted for the ground state according to the relation [jig 



a = (3 cot ^— (3^ 



which implicitly determines /3 as a function of the coupling constant a. 
In particular, this relation tells us that the critical coupling constant, 
ttc, must be approached in the limit (3^0. From the behaviour of the 
matrix elements {il)i\T\il)j) of the kinetic energy T in Eq. for large 
momenta p and of the matrix elements {il)i\Vc{r)\il)j) of the Coulomb 
interaction potential Vc{r) in Eq. (@) at small distances r, respectively, 
it should become rather evident that, for our choice of basis states iV^fc), 
these singularities will arise merely in those matrix elements which are 
taken with respect to l^/^o)? that is, only in (-i/^ol^l'^o) Q^nd ('0o|K:(''^)|'^o)- 
The remainder of our way is straightforward to go. With the above 
basis functions, we obtain for the matrix elements of the kinetic energy 
T in Eq. (p 

{^^\T\^3) 

2*+-?+2/5+i m T{i + (3+l) r(j + 13+1) 



X j dy 





TT ^T{2i + 2[3 + l)T{2j + 2(3 + l) 
cos[{i — j) arctan?/] — cos[(z + j + 2 P + 2) arctan?/] 



(1 + y2Yi+j+2p+l)/2 

which may be evaluated with the help of the expansion 

1 ^ f N\ fn7T\ 
cos(iVarctan^) = ^j^^ E ^ ^ j cos y 

for iV = 0,1,2,... 
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for the matrix elements of the Coulomb interaction potential Vc(r) in 
Eq. (0) 

^ ' ^ Vr(2^ + 2/? + i)r(2j + 2/? + i) 

and, finally, for the projections of the basis states \iljk) onto each other 

r(z + j + 2/3 + i) 



r(2^ + 2/3 + i)r(2j + 2/3 + r 



Of course, some care has to be taken when extracting the singularity of 
the matrix element (?/^o|7^|'^o)? which may be done by observing that 



lim 

/3^0 



, 1 - COS (2 + 2/3) arctand /• , 

ay , „ = 2 hm / ay 

y (1 + 7/2)1/2+/? /3^oy ^ 



2/ 



Introducing for notational simplicity a dimensionless and scaled energy 
eigenvalue e by 

^ 2 

-t/ = — m e , 

TT 

the characteristic equation for the Hamiltonian H in (0) thus becomes 





1 

41n2 - 


2 -£ 


V2 
3 




\ 




det 






17 


- e 




= 




3 




15 " 







V 



The roots of this characteristic equation are then given, for d = 1, by 

E = 2(21n2 - 1) , 
which entails the upper bound 

£"0 



= 0.4918. 



m 



and, for d = 2, by 



£ = -3- ("eO In 2 - 23 ± v/(601n2)2 _ 48001n2 + 1649) , 

J- ^ ^ 

which entails the upper bound 

^0 

— = 0.484288 . . . , 
m 
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for the ground-state eigenvalue £"0 of the Coufombic Hamiltonian 
In principle, the d (real) roots of the characteristic equation (g) may be 
determined algebraically up to and including the case d = 4, entailing, 
of course, analytic expressions of rather rapidly increasing complexity. 
For d = 4, this yields for the lowest-lying energy level the upper bound 

— = 0.4842564... . (5) 

m 

With respect to existing numerical determinations of energy levels, 
the exploration of the Hamiltonian (0) culminated so far with Ref. [jTl[ 



where — with the aid of the so-called local-energy theorem for the lower 
bound and a standard variational procedure for the upper bound — the 
admissible range of the energy eigenvalue of the ground state has been 
narrowed down to an amazingly tiny gap. In particular, at the critical 
coupling constant o^c this range is numerically fixed to be given by []TT 



^0 

0.4825 < — < 0.4842910 for a = . 

m 

Comparing these numerically determined bounds on the ground-state 
energy level with our above upper bounds resulting algebraically from 
the characteristic equation (|4|), one immediately realizes that, already 
for the case d = 2, our analytical bound lies well within the numerically 
obtained range. Consequently, the bound derived, in particular, in the 
case d = 4, Eq. (|5|), represents, rather obviously, a clear improvement 
of the best upper bounds known until now for the ground-state energy 
eigenvalue of the spinless relativistic Coulomb problem at the critical 
coupling constant. Needless to say, this simple analytical investigation 
serves, in addition, to strengthen our confidence in merely numerically 
computed upper bounds. 
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